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We use holography to study the spontaneous condensation of a neutral order parameter in a 
(2+l)-dimensional field theory at zero-temperature and finite density, dual to the electron star 
background of Hartnoll and Tavanfar. An appealing feature of this field theory is the emergence 
of an IR Lifshitz fixed-point with a finite dynamical critical exponent z, which is due to the strong 
interaction between critical bosonic degrees of freedom and a finite density of fermions (metallic 
quantum criticality) . We show that under some circumstances the electron star background develops 
a neutral scalar hair whose holographic interpretation is that the boundary field theory undergoes 
a quantum phase transition, with a Berezinski-Kosterlitz-Thouless character, to a phase with a 
neutral order parameter. Including the backreaction of the bulk neutral scalar on the background, 
we argue that the two phases across the quantum critical point have different z, a novelty that 
exists in certain quantum phase transitions in condensed matter systems. We also analyze the 
system at finite temperature and find that the phase transition becomes, as expected, second-order. 
Embedding the neutral scalar into a higher form, a variety of interesting phases could potentially be 
realized for the boundary field theory. Examples which are of particular interest to condensed matter 
physics include an antiferromagnetic phase where a vector condenses and break the spin symmetry, a 
quadrupole nematic phase which involves the condensation of a symmetric traccless tensor breaking 
rotational symmetry, or different phases of a system with competing order parameters. 



I. INTRODUCTION 



Condensed matter physics is replete with examples of 
ordered states described by a neutral operator acquir- 
ing a non-zero vacuum expectation value. Systems ex- 
hibiting magnetism [1], itinerant or otherwise, or more 
complicated ordered phases such as quadrupolar nematic 
states [2] are just a few examples of such ordered states. 
Our focus in this paper is to use holography to describe 
the onset of such neutral ordered states in itinerant-like 
d = 2 + 1 dimensional strongly-correlated systems with 
a finite dynamical exponent z. The standard method for 
treating the onset of magnetism in a Fermi liquid is to 
introduce a neutral bosonic order parameter at the level 
of a Hubbard-Stratonovich field in some fermionic model 
and then integrate out the fermions [3, 4]. Because the 
fermions are gapless, and hence inherently part of the 
low-energy degrees of freedom, they should not a priori 
be integrated out in the Wilsonian sense. Indeed, it is 
now well known [5-7] that this procedure is highly prob- 
lematic in two spatial dimensions. For example, in the 
case of d = 2 + 1 and z = 2, the fermions introduce an 
infinite number [5] of non-local marginal perturbations 
rendering any truncation of the resultant bosonic action 
at the Gaussian level moot. In the rcnormalization group 
sense, the critical theory in terms of the neutral bosonic 
order parameter does not possess any coupling constant 
from which a controlled expansion can be obtained. It 
turns out that even a 1/N expansion is of no use [6, 7] 
as the naive power counting fails, and an infinite number 
of planar diagrams contributes to the fermionic self en- 
ergy. Recent attempts [8] have been made to shed light 
on this problem by using holography where a strongly 
coupled quantum theory is mapped onto a weakly inter- 



acting dual theory of gravity . The boundary theory con- 
sidered in [8] is dual to the Reissner-Nordstrom AdS^+i 
black hole background (RN-AdS^+i) and flows in the IR 
to a (0+l)-dimcnsional CFT with a dynamical exponent 
z = 00. Hence, such a theory may have less direct bear- 
ing on real condensed matter systems, compared to a 
theory which flows in the IR to a fixed-point with finite 
z. Another drawback of this theory is the existence of a 
finite entropy density at zero temperature. 

All is not lost with holography in this context, however. 
Indeed, holography is ideally suited to treating this prob- 
lem as the failures of the existing conventional treatments 
[5-7] stem from the inherent strongly-coupled nature of 
the underlying theory. We show here that it is possible to 
construct a holographic setup in which a neutral bosonic 
operator condenses in an itinerant-like electron fluid with 
a finite value of z, and in d = 2 + 1. Interestingly, the case 
we consider also includes z = 2. A finite value of z in the 
IR is achieved by replacing the AdS 2 x near-horizon 
geometry of the RN-AdSd+i background with one that is 
of, say, Lifshitz type [13] as in [14]. Our work here then 
represents a holographic description of quantum critical- 
ity in a metallic system in two spatial dimensions with 
z = 2, and with a neutral bosonic order parameter. 

It was realized in [15] that a neutral scalar field with a 
mass square satisfying the Brcitenlohncr-Frecdman (BF) 
bound [16] of the asymptotic AdS^+i can cause an in- 
stability in the extremal, or near-extremal, RN-AdS^+i 
background if its effective mass square in the AdS 2 re- 
gion of the near horizon geometry violates the BF bound 



1 See [9—12] for reviews of holography with applications to 
strongly-correlated condensed matter-like phenomena. 
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of AdS2- As a result, the RN-AdS^+i background could 
become unstable to forming a neutral scalar hair 2 . The 
holographic interpretation of this hairy background is 
that the boundary theory is in a phase where the op- 
erator dual to the bulk neutral scalar field condenses. 
Such a condensate was explicitly constructed in [8] for 
the boundary theory dual to the RN-AdS^+i background, 
where it was shown that the operator condenses below 
some critical temperature as long as its conformal dimen- 
sion is in a certain range. The phase transition at finite 
temperature was shown to be second-order with mean- 
field exponents. Moreover, in the standard quantization, 
the critical temperature goes to zero as the conformal di- 
mension of the operator approaches a critical dimension 
from below. At zero temperature, on the other hand, 
by varying the conformal dimension of the operator, the 
boundary theory undergoes a quantum phase transition 
of the Berezinski-Kosterlitz-Thouless (BKT) type near 
the quantum critical point. The authors of [8] then used 
the condensation of the neutral scalar operator as a start- 
ing point to construct an antiferromagnetic phase in the 
boundary theory. This was accomplished by embedding 
the bulk neutral scalar field into a triplet whose dual op- 
erator condenses in such a way as to break the SU(2) 
spin symmetry to U(l) (see also [24]). The fluctuations 
of this antiferromagnetic order parameter (two Goldstone 
modes) were then analyzed in [8] and their dispersion re- 
lations were determined. 

As mentioned above, the boundary theory considered 
in [8] flows in the IR to a (0+l)-dimensional CFT with 
a dynamical exponent z — oo. This is an unwelcome 
feature in modeling metallic quantum criticality using 
holography. In this paper, we ask ourselves whether it is 
possible to condense a neutral scalar operator in a bound- 
ary field theory which flows in the IR to a fixed-point 
with finite z. We find that the answer is yes. To achieve 
this goal, one first replaces the AdS2 x K d_1 near horizon 
geometry of the RN-AdS^+i background with, say, a Lif- 
shitz geometry. Consequently, the background is asymp- 
totically AdSd+i while in the IR it has Lifshitz scaling. 
For simplicity, and with an eye towards condensed mat- 
ter applications, we assume that the boundary theory is 
d = 2 + 1 dimensional. A background with such proper- 
ties was recently constructed by Hartnoll and Tavanfar in 
[14] following [25, 26]. As in [14], we refer to this back- 
ground as the electron star background, or simply the 
star background. Compared to other backgrounds which 
are asymptotically AdS4 and Lifshitz in the IR [27], the 
electron star background has the appealing feature that 



2 This mechanism is identical to the one used for the construction 
of holographic superconductors [15, 17—19]. As discussed in [20, 
21], there is another way of making the AdS2 region unstable 
even if the effective mass square of the scalar field, either charged 
or neutral, satisfies the BF bound of AdS2- This mechanism 
involves deforming the boundary theory by relevant multi-trace 
operators [22, 23], which are, in this context, constructed out of 
the operator dual to the scalar field. 
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FIG. 1. A cartoon of the phase transition in our hologrphic setup 
as a function of the mass square, m 2 , of the bulk neutral scalar field 
4>. Once m 2 violates the BF bound, m 2 , of the far-interior Lifshitz 
region of the electron star background, the boundary theory goes 
into a phase where the dual operator <E> spontaneously condenses. 
The backreaction of the scalar on the geometry changes the value 
of the dynamical exponent z. The condensation of <t> is controlled 
by the IR characteristics of the background which remains d = 2 + 1 
Lifshitz across the transition. 



the IR Lifshitz scaling is due to the strong interaction 
between critical bosonic degrees of freedom and a finite 
density of fermions, which is of utmost importance in 
metallic quantum criticality. 

The (2+l)-dimensional boundary theory dual to the 
electron star background is a zero-temperature, but finite 
density, field theory which flows in the IR to a (2+1)- 
dimensional Lifshitz fixed-point with a finite z. We con- 
struct a neutral scalar hair for the electron star back- 
ground and show that upon varying the UV conformal 
dimension of the dual scalar operator, or any other knob 
in the boundary theory which allows the effective mass 
square of the bulk scalar in the far-interior Lifshitz region 
to change, the boundary field theory undergoes a quan- 
tum phase transition to a phase with a neutral order 
parameter. Figure 1 shows a cartoon of the phase transi- 
tion. A key characteristic we find is that the dynamical 
exponent changes across the transition. We analyze the 
nature of the phase transition at zero temperature and 
argue that it has a BKT character. We then consider the 
setup at finite temperature and show that below some 
critical temperature, the neutral order parameter con- 
denses and moreover the transition is of the standard 
mean field type. While the focus of this paper is mainly 
on the condensation of a single neutral scalar operator in 
the boundary theory, we discuss, very briefly, its appli- 
cations for the holographic modeling of phases with vec- 
tor or tensor order parameters, or systems with multiple 
order parameters. For example, following [8], an anti- 
ferromagnetic phase of the boundary theory in our case 
can easily be constructed by embedding the bulk neutral 
scalar into a triplet which is charged under the SU(2) spin 
symmetry, and allowing the dual neutral vector operator 
condense and break SU(2) down to U(l). Constructing a 
quadrupole nematic phase, where the order parameter is 
a symmetric traceless tensor which breaks the SO(2) ro- 
tational symmetry, is more difficult. We comment on the 
challenges of constructing such a phase using the embed- 
ding of the bulk neutral scalar into a symmetric traceless 
tensor. 

The paper is organized as follows. We start the next 
section by reviewing the electron star background of [14], 
and its finite-temperature variant [28, 29]. In section 
III we consider a neutral scalar field as a probe in the 
electron star background and show that under some cir- 
cumstances this gravitational system undergoes a phase 
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transition to a background with a neutral scalar hair. 
The holographic interpretation of the transition is that 
the boundary field theory undergoes a quantum phase 
transition of the BKT type from a Z 2 -symmetric phase 
to a phase with a neutral order parameter, which spon- 
taneously breaks this symmetry. We then consider the 
backreaction of the neutral scalar on the star background 
and show that the aforementioned two phases have differ- 
ent dynamical exponents. Section IV extends the analy- 
sis to the finite temperature case. In the last section we 
discuss various applications of our holographic setup and 
conclude with open questions and directions for future 
work. 



II. THE BACKGROUND 

The boundary theory that we wish to study is dual 
to the so-called electron star background. The zero- 
temperature background was constructed in [14], and was 
later extended to finite temperature in [28, 29]. In this 
section we briefly review the star background both at 
zero and finite temperature, and set the stage for later 
calculations. 



A. Zero- Temperature Background 

Consider an ideal fluid of charged fermions in a (3+1)- 
dimensional space-time with a negative cosmological con- 
stant A = —3/L 2 . The equations of motion of the system 
read 

Rmn ~ -^RguN + A^mjv = k 2 Tmn, (1) 
V M F MJV = e 2 J N , (2) 
where n 2 = 8nG^, e is the fermion charge and 

F PQ F p Q) 



Tmn — ~2 {FmpF^ — -gMN 

+ (p + p) u M u N +pg M N, 
Jm = & um, 



(3) 
(4) 



with p, p, a and u a being, respectively, energy density, 
pressure, charge density and four-velocity (satisfying the 
normalization umu m = —1) of the fluid. Throughout 
the paper, we denote the bulk indices by capital letters 
M, N, ■ ■ ■ = {t, x, y, r} while the boundary theory direc- 
tions are denoted by the Greek indices p, u, . . . = {t, x, y}. 

To solve the above equations, we take the ansatz for 
the metric and the gauge field to be of the following form 



ds 2 
L 2 



-f(r)dt 2 



1 



eL 

A = — h(r)dt, 



(dx 2 +dy 2 ) +g(r)dr 2 , (5) 

(6) 



where the asymptotic boundary is at r — > 0. In ad- 
dition, an equation of state for the fluid must also be 



specified. Suppose the fluid is composed of zero temper- 
ature charged fermions with mass mf, with the density 
of states 



n(E) = f3EJE 2 - m 2 



(7) 



where f3 is a constant of order one. The authors of [14] 
then assumed a locally flat approximation where, at each 
r-slice, the fermion dynamics are determined by the local 
chemical potential 



e h(r) 

*7wy 



(8) 



and argued that such an approximation is self-consistent 
in a regime of parameters where gravity can be treated 
classically with an order-one backreaction of the fermions 
on the geometry. Given this approximation one can 
then use (7) to compute the fluid energy density, pres- 
sure and charge density as function of r. Obviously, 
n(E) = p(r) = p(r) = a(r) = for r < r s where r s is the 
boundary of the fluid given by solving p\ oc (r s ) = mf. 

It is more convenient to rescale p —> p/(L 2 K 2 ), p — > 
p/(L 2 n 2 ), a -> cr/(L 2 en), f3 -> ( K 2 /e 4 L 2 )p, m f -> 
(e//t)mf and p\ oc —> (e/fc) pi oc , with the new quantities 
being all dimensionless. Putting everything together, the 
equations of motion then take the form 

1 (fir) g'(r) 4\ g(r)h(r) . . 

+ + - ) + /^X' a(r) = 0, (9) 



r \ f( r ) g( r ) r 

f(r) h'{rf 



+ [3 + p(r)]g(r) 



0, (10) 



r/(r) 2/(r) 
h"(r) + (Ih(r)h'(r) /(r)) a(r) = 0. (11) 

In the far-interior region r — > oo, the solution is of the 
Lifshitz type [14]. As argued there, in order for the Lif- 
shitz solution to make sense, one should have < mf < 1. 
The solution in the interior region r s < r < oo (where the 
fluid energy density, pressure and charge density are all 
non-zero) can easily be worked out by perturbing the Lif- 
shitz solution, and demanding the perturbation to grow 
as r — > 0, but not blow up as r — > oo. In so doing, one 
finds 



1 



f(r) = ^(l + kr a - +■■■), 
g(r) = ^(l + Ql r a - +•••), 
h(r) = 1(1 + ^" +■■■), 



(12) 
(13) 
(14) 



where z = z(m,f, (i) is the dynamical critical exponent 



2 = 



36z 4 (z - 1) 



[(l-m 2 )z-ir 



(15) 



and 0i, \)i and all higher order perturbation coefficients 
are determined in terms of f i , which itself could be set to 
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any value by rescaling r, t, x, y. Also, a_ is determined in 
terms of m.f and z, the expression of which can be found 
in [14] . The dependence of z on mf and f5 is complicated, 
but satisfies z> (1 — m 2 ) -1 > 1. (A numerical plot of z 
in terms of /3 for some sample values of mf, along with 
the explicit forms of its asymptotic expressions can also 
be found in [14].) The solution in the exterior region < 
r < r s (with the fluid energy density, pressure and charge 
density all being zero) is RN-AdS4 where the metric and 
gauge field functions are given by 

/W = ^(c 2 -Mr 3 + iQV), 9{r) = ^r y 

h(r) =fi-Qr. (16) 

Matching /(r), g(r), h(r) and h'(r) of the interior and 
exterior solutions at r = r s will then determine the con- 
stants c, M, Q and /i. Note that the constant c is intro- 
duced because choosing a particular value of fi in the in- 
terior solution fixes the normalization of time. The elec- 
tron star background is a family of asymptotically AdS4 
solutions parametrized by mf and fj. The far-interior re- 
gion of the background is generically Lifshitz. But, in 
some limits of the fj-m,f parameter space, the far-interior 
region is characterized by a geometry other than Lifshitz. 
For example, an AdS2 x R 2 solution (where z = oo) can 
be recovered either by taking /3 — > while keeping mf 
fixed, or by taking the mf -> 1 limit from below. An- 
other interesting limit of the far-interior solution is AdS4 
(with z = 1) which is obtained by taking j3 — > oo at 
mf = 0. 

In this paper, we exclude the two limiting AdS2 x R 2 
and AdS4 interior solutions by taking mf e (0,1). In 
other words, the far-interior solution is Lifshitz with a 
finite z > 1. In this case, one obtains q > 1, which turns 
out to be crucial in our subsequent discussions. Note that 
the zero-temperature electron star background has zero 
entropy density. 



B. Finite- Temperature Background 



the other hand, if mf is small enough, there exists a crit- 
ical radius r = r c for which mf equals the local chemical 
potential /j\ oc at that radius. This happens at a critical 
temperature T = T*. The quantities r c and T* arc de- 
termined by solving the two equations y(ii oc (r s ) = mf and 
Moc( r s) = 0> w hh prime denoting derivative with respect 
to r. 

For T < T*, there exists a finite width r 2 > r > n 
for which the fluid has non-zero energy density, pressure 
and charge density. The radii r\ and r 2 are solutions to 
Mioc(?"i,2) = mi. Thus, the finite-temperature electron 
star background is divided into three regions. In the so- 
called inner region of the background, r > r > r 2 , the 
solution is a RN-AdS4 black hole with 

h(r) = Mo - Qor, (17) 

where ro is the horizon radius given by the largest real 
root of f(r) = 0, Mo = Qo^o and M r% = 1 + Mo r o/ 2 - 
The temperature of the black hole is given by T = 
(47tc) _1 \df(ro)/dr\, where the parameter c is introduced 
for later convenience. The quantity Mo is the inner region 
chemical potential which is proportional to the charge Qq 
hidden "inside" the black hole horizon. Note that n does 
not correspond to the chemical potential of the bound- 
ary field theory. Indeed, can be used as a tuning 
parameter to change the temperature and the chemical 
potential of the boundary field theory. In the intermedi- 
ate region r 2 > r > r\, the fluid pressure, energy density 
and charge density are all non-zero. The background in 
this intermediate region is obtained by solving the equa- 
tions (9)— (11). In the exterior region, namely r x > r, the 
background is again described by a RN-AdS4 solution 
with 

/W = ^(l-Mr 3 + iQV), ff(r) = ^, 
h(r)=c( f i-Qr), (18) 



The details of the finite-temperature version of the 
electron star background are as follows [28, 29]. Tem- 
perature is introduced by making the Euclidean time 
coordinate periodic with a period 1/T, where T is the 
temperature. This then requires a horizon at some finite 
radius r = tq in the interior of the geometry. As argued 
in [28, 29], in the regime of validity of the electron star 
background e <~ k/L <C 1, the effect of temperature on 
the fermion density of state could be neglected. There- 
fore, the density of states is given by (7) and with the 
same planar ansatz for the metric and the gauge field as 
before, one has to solve the equations (9)— (11) in order 
to determine the finite temperature background. Again, 
for large enough mf, the fluid energy density, pressure 
and charge density vanish everywhere. As a result, a 
RN-AdS4 background will satisfy those equations. On 



where the boundary theory quantities c, M, Q and /i are 
determined by matching the solutions at r = r\ and r 2 . 
Note that \x is the chemical potential of the boundary 
field theory. 

In summary, for mf € (0, 1), as the temperature is low- 
ered, the RN-AdS4 background undergoes a phase transi- 
tion to a finite-temperature electron star background. It 
was shown in [28, 29] that this phase transition is third or- 
der. Also, it is worth emphasizing that an appealing fea- 
ture of the finite-temperature electron star background 
is that, unlike the RN-AdS4 solution, the entropy den- 
sity s approaches zero as s ~ T 2 / z in the T -> limit 
[29], where, for given allowed values of the parameters 
mf and /?, the dynamical critical exponent z agrees with 
the one obtained from the zero-temperature electron star 
background with the same values of mf and (3. 
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III. NEUTRAL SCALAR ORDER PARAMETER 
AT ZERO TEMPERATURE 

In this section we analyze the spontaneous condensa- 
tion of a neutral scalar operator in a strongly-coupled 
(2+l)-dimensional zero-temperature boundary field the- 
ory with a finite density of fermions, which is dual to the 
the zero-temperature electron star background reviewed 
in the previous section. In the context of the AdS/CFT 
correspondence, the question of whether a scalar oper- 
ator can spontaneously condense in a strongly-coupled 
boundary theory maps to the question of whether the 
gravitational background can develop a scalar hair. So, 
we are naturally led to investigate whether the electron 
star background can develop a neutral scalar hair. 



The action (19) is invariant under a Z2 symmetry: 
— > —(f). This symmetry is broken in a phase of the 
boundary theory where the operator dual to (f> sponta- 
neously condenses. The symmetry group can easily be 
enlarged by considering more scalar fields, and depend- 
ing on which dual operators condense, a variety of inter- 
esting symmetry breaking patterns could be realized in 
the boundary theory. Such holographic setups can po- 
tentially be used to capture the relevant physics of the 
strongly-correlated systems with competing order param- 
eters. In this section, we deal with just one neutral scalar 
field with an action given in (19). 

The equation of motion for <j> reads 



dm (g MN V=gd N <f>) = V'{<f>), 



(21) 



A. Neutral Scalar Hair for the Electron Star 

We introduce the scalar field <j> in the electron star 
background by considering the action 



2k 2 X 



1 



(V0) 2 + V(<t>) 



where the potential V(<j>) takes the form 



V(<P) 



1 

4L2 



+ m 2 L 2 y-m 4 L 4 



(19) 



(20) 



We choose the mass square of the scalar, m 2 , to be neg- 
ative but above the BF bound of the asymptotic AdS 4 
region. Namely, the mass of the scalar is assumed to be 
in the range —9/4 < m 2 L 2 < 0. Note that since m 2 
is taken to be negative, the field <f> should be expanded 
asymptotically around the maximum 3 of the Mexican- 
hat potential V{4>), i-e. <f> = 0. As we will explain below, 
it is important that V(<fi) has (at least) a minimum [8]. 
Note that the constant term in V((f>) is there to ensure 
zero contribution to the cosmological constant of the elec- 
tron star background, A = — 3/L 2 , for (f> = 0. Also, A is 
a coupling constant chosen, at first, to be large in order 
for the probe approximation to be valid. We later re- 
lax this assumption when wc consider the backreaction 
of the scalar on the electron star background. The same 
scalar action was used in [8] to study the condensation of 
a neutral scalar order parameter in the boundary theory 
dual to the RN-AdS4 background, both at zero and finite 
temperature 4 . 



3 Expanding <j> in the asymptotic region around either of the two 
minima of V(tj>) corresponds to a positive m 2 . We will not con- 
sider such a case in this paper, as, with the mechanism explained 
in the paper, we do not expect to see a condensation of the dual 
neutral scalar operator in the boundary theory. 

4 See also [21] for a discussion of condensing neutral scalar opera- 
tors in holographic setups using multi-trace deformations. 



where prime here denotes derivative with respect to (ft. 
Given the background metric in (5) and choosing the 
ansatz (ft — <ft(r), equation (21) takes the form 



r Q M 




L 2 V'{(j>). (22) 



We are after a non-trivial solution which is regular as 
r — > 00 and normalizable as r — »■ 0. Note that r — > 00 
is the far-interior of the background which is Lifshitz, 
while r — > is the boundary of the background which 
is AdS4. From (12) and (13), the metric functions in 
the interior region take the form f(r) = r~ 2z + ■ ■ ■ and 
g(r) = Qr~ 2 + • • • , where the dots represent the sub- 
leading terms as r — > 00. Hence, the equation of motion 
for (ft in the far-interior region reduces to 



n 3+z 







8 r [r-^d r ^j = L 2 V\4>). 



(23) 



Demanding a non-trivial regular solution in the far inte- 
rior forces (ft to sit at a minimum of the potential V{<ft). 
This can easily be shown following [8] . Suppose the non- 
trivial regular solution in the far interior of the back- 
ground behaves as (ft — <fto + S(f> where 8<ft ~ r s with 8 < 0, 
so that (f> — > (f> as r — > 00. One then easily deduces that 
the left-hand side of the equation (23) vanishes, resulting 
in V'((f)o) = 0. So, the scalar field must sit at an ex- 
tremum of the potential. Now linearizing equation (23) 
around (ft = (f> n , we find 

-r z+3 d r (r- {1+ ^d r 8<p) - L 2 V" {(}> )8cf> = 0(8(f> 2 ). (24) 

Q V J 

The general solution of the above equation takes the form 
8(f){ r ) = ar s + + br s - with 



S± 



(2 + z) ± v/(2 + z) 2 +4 W(</> o ) 



(25) 



Because the solution ar s + blows up in the r — > 00 limit, 
we set a = 0. In order for br 5 ~ to be a regular solution 
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FIG. 2. Plots of B versus A for the neutral scalar field <j> with a mass square m 2 L 2 = —2.20 on an electron star background with 
nif = 0.36 and (a) z = 2 and (b) 2 = 3. The BF bounds in the far-interior of the star (with mf = 0.36) for 2 = 2 and 2 = 3 are 
m 2 L 2 ~ —2.12 and —2.02. respectively. Since B is non-zero at A = 0, the plots show the spontaneous condensation of the dual operator 
<3? in the respective boundary theories. 



in the interior, one has to have V"(<fio) > (note that g 
is positive). This then proves that the scalar field must 
sit at a minimum of V((j>), if one demands a non-trivial 
regular solution for cj>(r) as r — > oo. Thus, the regular 
solution takes the following form in the far-interior region 
of the background 



(r) = 0o + br s - (! + •••) 



r — > oo, 



(26) 



with 



5- = 



(2 + z)- V(2 + ^) 2 - 



i 2 L 2 



(27) 



Asymptotically as r — ► 0, from (16), we have f(r) 



' + ■ ■ ■ and g(r) = r + • 



The scalar equation of 



motion (21) then becomes 
1 



+ m z L z ) = 0, r^0. (28) 



Linearizing equation (28) around = 0, the general so- 
lution will be of the form 



^ = ir J - A (l + ---) + Br a (l + .-.), 



(29) 



where A = 



3 



rn 



2 L 2 



, v ■ ■ !■ For -| < m 2 L 2 < -f , 
either A or £? could be called the source and the other 
one the vacuum expectation value (vev) of the dual op- 
erator $ [30]. We choose to work in the so-called "stan- 
dard quantization" where we take A to be the source for 
<I>, and B as the vev of $, denoted by ($). (Note that 
for m 2 L 2 > — |, only the standard quantization exists.) 
Hence, A is the dimension of the neutral scalar opera- 
tor $, which is the holographic dual of the bulk neutral 
scalar field (j>. Given the boundary condition (26) for <fi in 
the interior, for each value of b, the equation (22) could 
be numerically integrated to the asymptotic r — > region 
from which both A and B could be determined. Repeat- 
ing the process for different values of b, the dependence of 
A and B on b, and thus the dependence of B on A, could 




FIG. 3. Plots of m\l? = -(2 + 2) 2 /4g as function of z for 
mf = 0.36 (red) and mf = 0.7 (blue). Here L is the curvature 
radius of the asymptotic AdS4 region. As z increases the curves 
asymptote to the dashed line, which represents the scalar BF bound 
in AdS 2 . 



easily be obtained. Since we are interested in the con- 
densation of the operator $ without turning on a source, 
we look for solutions where B remains finite at A = 0. 

In performing the numerics in this section, it is more 
convenient to define a dimensionless radial coordinate by 
rescaling r — > r s r. As a result, we also make 6, A and 
B dimensionless by redefining b — > r s S ~ b, A — > r^~ 3 A 
and B — > r~ A B. We also set L = 1 in our numerics 
throughout this paper. Figure 2 shows B versus A for 
the bulk neutral scalar with m 2 L 2 = —2.2 on two dif- 
ferent electron star backgrounds, one which flows in the 
IR to a theory with z = 2 and the other to a theory 
with z = 3. For both of the star backgrounds, we set 
rnf = 0.36. As the plots in Figure 2 demonstrate B is 
non zero when A is zero, indicating that the neutral op- 
erator $ spontaneously condenses in the two boundary 
theories (which are the holographic duals of the two elec- 
tron star backgrounds). Said another way, the aforemen- 
tioned two electron star backgrounds have each formed a 
neutral scalar hair. 

The reasoning behind the formation of a neutral scalar 
hair for the electron star background is similar to the RN- 
AdS4 case discussed in [8, 15]. Recall that the (2+1)- 
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dimensional boundary field theory under consideration 
flows in the IR to a (2+l)-dimcnsional Lifshitz fixed- 
point with a finite dynamical critical exponent z. Let's 
denote the dimension of the operator $ at the IR Lifshitz 
fixed-point by Air. One can then easily show that 



A IR = ^ 



(2 + z) + v /(2 + z) 2 +4 m 2 L 2 



(30) 



In the discussion below equation (20), we have assumed 
that the mass square of <f> satisfies — | < m 2 L 2 < 0, where 
the lower bound is the BF bound of the asymptotic AdS4 
region of the electron star background. Since for the IR 
Lifshitz fixed-point with z > 1 one has g > 1, there is 
always a range of mass square for which the expression 
under the square root in (30) becomes complex, while still 
satisfying the BF bound of the asymptotic AdS4 region. 
Namely, for 



< m 2 L 2 



< m 2 L 2 , 



with m n L 



1 T 2 



4g 



(31) 



the bulk neutral scalar field becomes tachyonic in the far- 
interior region or equivalently, in the language of the the 
dual field theory, the scaling dimension of $ at the IR 
Lifshitz fixed-point Air becomes imaginary, signaling a 
potential instability in the theory. Note that — (z+2) 2 /4g 
is the scalar BF bound in the far-interior Lifshitz geom- 
etry of the electron star background. Figure 3 shows 
plots of this bound as a function of z for mi — 0.36 
(red) and mf = 0.7 (blue). One sees from the plots 
in Figure 3 that as z increases, m 2 L 2 approaches the 
value -3/2, which equals the scalar BF bound in an AdS2 
background. (Note that L is the curvature radius of the 
asymptotic AdS4 region.) Such behavior is expected, and 
is due to the fact that for a fixed value of mf € (0, 1) and 
in the z — > oo limit, the far-interior region of the elec- 
tron star background becomes AdS2 x IR 2 [14]. In other 
words, the electron star background reverts back to the 
RN-AdS 4 solution. 

In both of the plots in Figure 2, the mass square of the 
bulk scalar m 2 L 2 = —2.2 violates the BF bound of the 
far-interior Lifshitz region of each of the two star back- 
grounds while satisfying the BF bound of their asymp- 
totic AdS 4 regions. Hence, we see from the plots that 
the neutral scalar operator $ spontaneously condenses 
in each of the two boundary theories, and breaks the 
aforementioned Z 2 symmetry. On the other hand, if the 
mass square is above the BF bound in the far-interior 
Lifshitz region, no spontaneous condensation for $ is 
seen 5 in the boundary theory, and, as a result, the Z 2 
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FIG. 4. B versus m 2 L 2 for the neutral scalar field <j> on an 
electron star background with mf = 0.36 and 2 = 2. Due to 
numerical limitations we could not access the region of mass square 
very close to the BF bound, m 2 L 2 ~ —2.12, of the far-interior of 
the star. The inset shows a close up of the behavior of B for larger 
values of m 2 L 2 . 



symmetry is intact. The plot in Figure 4 shows the be- 
havior of B as a function of m 2 L 2 for the bulk neutral 
scalar in an electron star background with m.f = 0.36 
and z = 2. As the plot shows, B approaches zero as 
m 2 L 2 approaches m 2 L 2 ~ —2.12 from below. In general, 
by tuning the mass square of the bulk scalar across the 
critical value m 2 L 2 , we obtain a hairy electron star back- 
ground. Or, simply put in the language of the boundary 
theory, varying the dimension A of the operator <E> across 
a critical dimension A c , the boundary theory at low en- 
ergies goes through a quantum phase transition, from 
an uncondensed (Z 2 -symmetric) phase, ($) = 0, to a 
phase with a neutral scalar order parameter, ($) 7^ 0, in 
which the Z 2 symmetry is spontaneously broken. Here, 

A c = i (3+v/9-(z + 2) 2 /a). 

A few comments are in order. We would like to em- 
phasize that the inequality q > 1 plays an important role 
here. Indeed, should q < 1, the BF bound in the interior 
Lifshitz region, for the case of z > 1, would be less than 
the scalar BF bound in the asymptotic AdS 4 region, and, 
as a result, there would be no possibility for the phase 
transition of the kind we just described above. 

The two phases described above have the same dy- 
namical critical exponent z. This is just an artifact of 
the probe approximation we have employed so far where 
the backreaction of the bulk neutral scalar field on the 
electron star background in neglected. In a subsection 
to follow, we will take into account such a backreac- 
tion and show that the aforementioned zero-temperature 
quantum phase transition is indeed between two phases 



5 Adding a double-trace deformation of the type g J d?x $ 2 (with 
g < 0) to the boundary theory action, and choosing the alterna- 
tive quantization where the dimension of the operator <I> satisfies 
1/2 < A < 3/2, makes it possible for <3? to condense even if 
the mass square of the dual bulk neutral scalar is above the BF 
bound of the far-interior Lifshitz region. It would be interesting 



to explore, along the lines of [21], the consequences of the afore- 
mentioned double-trace deformation for the holographic setup 
under consideration here. In the standard quantization, how- 
ever, the mass square of the neutral scalar being below the BF 
bound of the far-interior Lifshitz region is the necessary and suf- 
ficient condition for the condensation of <1> in the theory dual to 
the zero-temperature electron star background. 
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with different z. This is a particular novelty in certain 
quantum phase transitions in condensed matter systems. 
Phase transitions in which the critical dynamical expo- 
nent changes across the critical point have been suggested 
to underlie magnetism in the cupratcs [31, 32] and in 
nematic-smectic transitions [33]. 

The existence of the neutral scalar hair is essentially 
goverened by the physics in the far-interior region of the 
electron star background. More concretely, in the stan- 
dard quantization, it is the effective mass square of the 
neutral scalar field in the Lifshitz region which deter- 
mines whether the electron star background becomes un- 
stable to forming a neutral scalar hair. Since the bulk 
scalar field is neutral, its effective mass in the Lifshitz 
region is the same as the mass m in the asymptotic AdS 4 
region. So, tuning m 2 L 2 across the BF bound in the 
Lifshitz region, m 2 L 2 = —(z + 2) 2 /4g, the electron star 
background undergoes a phase transition to the forma- 
tion of a neutral scalar hair. Tuning the mass square of 
the scalar in the asymptotic region is not the only way 
to achieve this phase transition. For example, a simple 
way for tuning the effective mass square of the neutral 
scalar field in the far-interior Lifshitz region is to couple 
the scalar to the square of the Weyl tensor in the action 
(19), by adding to the potential V{4>) a term of, say, the 
following form [34] 



SV 



L 2 l 2 6 2 W 2 



(32) 



where W 2 = W M npqW mnpq with W M npq being the 
Weyl tensor. The effect of this new coupling is to make 
the mass of the scalar to depend on the radial coordinate 
r. Asymptotically, W 2 — > as r — > 0. So the mass of 
the scalar in the asymptotic AdS 4 region is untouched 
by the coupling (32). In the far-interior region, one can 
easily show that W 2 -> 4z 2 (— l + z) 2 /(3g 2 L 4 ) as r -> oo, 
hence, the mass square of the scalar in the Lifshitz region 
is shifted by a constant proportional to I 2 . Thus, the ef- 
fective mass square of the scalar in the Lifshitz region 
can be varied simply by varying I. The coupling I could 
be tuned such that the effective mass square in the Lif- 
shitz region violates the BF bound there, in which case 
it is not hard to show that the electron star background 
develops a neutral scalar hair. Neither tuning to nor tun- 
ing I is an operation that could be performed in a single 
theory. If one thinks of backgrounds such as the electron 
star background, or RN-AdS4, as coming from, say, M- 
theory, then the couplings to, /, • • • are generally fixed. 
In this regard, tuning m or I seems unnatural. Never- 
theless, we continue to tunc m 2 L 2 in order to be able 
to obtain the aforementioned phase transition. As the 
transition owns its existence to what the effective mass 
square of the scalar is in the far-interior region, we do 
expect the qualitative features of the phase transition to 
stay the same regardless of whether it has been obtained 
by tuning m, I, or a UV knob (such as a magnetic field) 
that could be continuously changed within a single the- 
ory. A simple way to introduce such a knob that can be 



continuously varied in a single theory has been explained 
in [8]. 



B. Free Energy of the Condensed Phase 

In this section we determine the free energy density of 
the boundary theory in the condensed phase ($) ^ 0. 
The free energy is given by the on-shcll action, up to 
boundary counter terms. Assuming the bulk scalar field 
does not back-react on the background, the free energy of 
the condensed phase is just the sum of two contributions: 
one coming from the electron star background and the 
other from the neutral scalar field. We work in the grand 
canonical ensemble where the boundary theory chemical 
potential \x is held fixed. As shown in [14], the free energy 
density of the theory dual to the zero-temperature elec- 
tron star background is given by f2i = M — /j,Q — —M/2, 
with M and Q being the energy and charge densities, re- 
spectively. Thus, we only need to determine the free 
energy density contributed by the boundary theory op- 
erator 3>, denoted hereafter by £1 2 - 

Up to boundary counter terms, f2 2 is simply given by 
the action (19) evaluated on-shell. To make the expres- 
sions less cluttered, we drop the coefficient 1/(2kA) from 
the action. We also set L = 1 in this discussion. Upon 
integrating by parts and using the equation of motion 
(21), the action (19) takes the form 



! d *-{ 



-99 



+ 



(33) 



where prime denotes derivative with respect to r and e is 
a cutoff introduced to regulate the above integral, as it 
is naively divergent. We will eventually send e to zero. 

If we work in an ensemble where the source is kept 
fixed, i.e. 8<f> — 0, the on-shell action (33) needs to be 
supplemented by a counter term 



;(3-A) 



(34) 



to cancel the divergence. Note that 7 is the induced 
metric on the boundary. Adding the above counter term 
to (33), and taking the e — > limit, the renormalized 
on-shell action S Ten becomes 



V 



-(2A- 3)AB + 



1 f 

4Jo 



dr\ 



(35) 



where V = J d 3 x, and we have substituted the asymp- 
totic expansion of the scalar field at the boundary. Since 
O2 = — Srcn/V is the free energy density contributed by 
the operator the free energy density of the system 
fi = r^i + 2 then reads 



fi = - 



M 



-|(2A 



1 r 



(36) 
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FIG. 5. Plot of A!)//i 3 versus m 2 L 2 . For the background, we 
set z = 2 and m,{ = 0.36. Numerical difficulties prevented us from 
obtaining more data points for the values of m 2 L 2 very close to 
the critical value m 2 L ~ —2.12. 



If instead wc work in an ensemble where the response 
is fixed, i.e. 5<j)' = 0, then one has to add an additional 
counter term, analogous to the Gibbons-Hawking term 
[35, 36], to the on-shell action in order to render the 
variational problem well-defined. The total counter term 
for this case becomes 



+ 



J d 3 x 
|(3-A)/ 



(37) 



Going through the same steps as before, we write the free 
energy density of the system as 



M c. nA 

T + 2 (2A 



Z)AB 



4 Jo 



(38) 



The condensed phase, ($) ^ 0, corresponds, in the 
standard quantization, to A = and fi / while the 
uncondensed phase is given by A = B = 0. Indeed, in 
the latter phase, there is only the trivial solution <j>(r) = 
in the bulk. The change in free energy is then given by 



n 



condensed 



uncondensed 



(39) 



Therefore, we see that when the operator $ condenses, 
the free energy of the system is actually lower than that 
of the uncondensed phase. Figure 5 shows a plot of 
the change in free energy density, AQ — ri C ondonscd — 
^uncondensed, versus m 2 L 2 . Clearly, the condensed phase 
has a lower free energy than the uncondensed phase, the 
difference approaches zero when the critical point is ap- 
proached from below. 



v is a parameter that can be tuned across the critical 
point and c\ > 0. 

Consider the boundary theory dual to the AdS^+i 
background, with a UV scale. For this theory, the au- 
thors of [37] have argued on general grounds that if the 
mass square of a bulk scalar is slightly below the BF 
bound, an IR scale A IR will be generated, which is related 
to the UV scale Auv via Air = Auv exp(— c-ij \Jv c — v) 
with C2 > 0. Concrete examples of this phenomenon were 
constructed in [8, 38] using a background with an AdS2 
near horizon geometry. Also, it was argued in [38] that 
the same phenomenon would occur in a boundary theory 
dual to a Lifshitz-like background [38] . As shown below, 
the discussion in [8, 37] for a BKT-type phase transition 
can handily be generalized to our case as well. 

Consider the equation of motion for <j> given by equa- 
tion (23). Linearizing around (f> — in the far-interior 
Lifshitz region of the electron star background, this equa- 
tion can be put in the form 



X"(r) 







-r 2z co 2 



+ ( 



TO 



ml)L 2 



1 

4fl 



X(r), (40) 



where we assumed a time-dependence of e~ tult for <p, and 
defined <p{r) = A z+1 ^ 2 x(r). Note that in the above equa- 
tion m 2 L 2 = -(2 + z) 2 /Aq is the BF bound of the far- 
interior Lifshitz region. For the mass square below the 
BF bound, there will exist a negative-energy bound state 
if < r < oo. On the other hand, if the geometry is cut 
by some suitably chosen UV and IR "walls" , denoted by 
ttjv and hr such that hr > r > ruv, the existence of 
negative-energy bound states could be prevented. This 
is actually the case in our setup. The far-interior Lif- 
shitz region does not extend all the way to the boundary, 
and at some radius matches to the exterior geometry. 
Also, the IR wall is provided by the condensate, whose 
existence is essentially due to the stabilizing effect of the 
quartic term in the scalar potential (20). 

Since we are interested in the onset of an instability 
we look at the lo = solutions. As argued in [8, 37] the 
exact form of the boundary conditions for <j> at the walls 
arc not important and could be taken to be Dirichlet. 
So, we choose </>(ruv) = 4>{ r m) = 0. The zero-energy 
solution satisfying the Dirichlet boundary condition at 
ruv is found to be [8, 37] 



cj)(r) = r (*+ 2 )/ 2 sin 



V Q (m 2 L 2 — m 2 L 2 ) log ■ 



■ (41) 



To satisfy the boundary condition at hr , one should then 
have 



Zero-Temerature Phase Transition 



rm = ruv exp 



\Jq (m 2 L 2 — m 2 L 2 ) 



(42) 



In this subsection we analyze the nature of this zero- 
temperature phase transition and argue that it has a 
BKT character, meaning that close to the critical point in 
the condensed phase ($) <~ [i A exp(— c\ j^Jv c — v) , where 



Note that HR/ruv = Atjv/Air. For hr given by the 
above formula, the ground state has zero energy. So, 
there will be no negative-energy bound states. On the 
other hand, for larger rm there will be a negative-energy 
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ground state, hence instability. For the mass square of 
the dual bulk scalar slightly less than the BF bound, the 
IR Lifshitz fixed-point is still scale invariant over a large 
energy scale. For energies below Air, the operator $ 
wants to condense and end the instability. For the mass 
square slightly less than the BF bound, the IR (mass) 
dimension of $ according to (30) is almost (z + 2)/2. 
So, following the argument given in [8], and given that 
Ayv ~ P, one then concludes 



<<*>> 



p A cxp 



(z + 2)tt 



2^0 (m 2 c L 2 -m 2 L 2 ) 



(43) 



Note that there are also an infinite number of states, 
characteristic of Efimov states [8, 38, 39], with (<&)„ ~ 

p, A exp — (1 + z/2)nir / yj q (m 2 L 2 — m 2 L 2 ) , where n = 

2, 3, 4, • • • . The ground state, however, is given by (43). 
The scaling of the condensate as in (43) signifies that the 
underlying transition is of the BKT-type. 



D. Backreaction 

In this section, we consider the backreaction of the 
neutral scalar field on the electron star background, and 
show that if the mass square of the scalar is below the BF 
bound of the far-interior region, the dual operator in the 
boundary theory condenses. Moreover, in the condensed 
phase the value of the dynamical critical exponent z is 
different from that of the uncondensed phase. 

To consider the effect of backreaction of the scalar field 
on the background, we need to modify the expression for 
Tmn in (3) by including the energy- momentum tensor 
of the scalar. Note that since the scalar field is neutral, 
the current Jjv is unchanged. The ansatz that we take 
for the metric and the gauge field is the same as in (5) 
and (6), respectively. For the scalar, as in the previous 
sections, we take the ansatz <j> = <j)(r). Substituting the 
ansatz into the Einstcin-Maxwell-fluid-scalar equations 
of motion, one obtains 



= 



KM + 7M + ?) + b(r) + * )]9(r, + ^' <r)2 ' 



1 f'(r) h'(r) 
= P'(r) + - ]p(r) + p(r)\ - (7 ( r)) 



f{r) h'{rf 



rf(r) 2/(r) 



[3+p(r)}g(r) 



4A 



4>'{r) 



1 9J £ L2V ^ 



= h"(r)+ r -[p(r) + p{r)\ g(r)h'(r)-g(r)y/f(r)v(r) + ^rti {r)<j>' {r)\ 



(44) 
(45) 

(46) 
(47) 
(48) 



As before, to solve the above equations, we need to 
specify the equation of state for the fermion fluid. Com- 
pared to the case without the scalar backreaction, the 
equation (45) is unchanged . This is due to the fact that 
in our model, the neutral scalar does not couple directly 
to the fermion fluid. Thus, the fermion density of states 
is given by (7). Also, similar to the case without the 
scalar backreaction, we adopt a locally flat approxima- 
tion where the fermion dynamics is determined by the 
local chemical potential, whose form is given in (8). As 
a result, the formulae for a, p, and p are the same as 
before, with the understanding that the metric and the 
gauge field functions in these expressions are now the 
backreacted ones. We do not write the formulae for a, p, 
and p here as they can be found in [14] . It can easily be 
verified that the expressions for cr, p, p satisfy the equa- 
tion (45). Thus, in order to determine the backreacted 
background, one has to solve the remaining equations, 
namely the equations (44), (46), (47) and (48). 

In the far- interior region r — > oo, we anticipate a 



Lifshitz-type geometry. This is due to the fact that the 
local charge density screens the electric field making it 
massive, and massive vector fields usually give rise to 
Lifshitz geometries [13]. Indeed, in the region r — > oo, 
the following 



L 2 



+ l(rf.T 2 +dj/ 2 ) + 4dr 2 , 



A = -dt, 

k r z 



(49) 



is an exact solution to the equations of motion, with (j> — 
0, ±\/—m 2 L 2 . The expressions for q and 1} are identical 
to the ones given in (15) except that, with the scalar 
backreaction included, and when (f) = ±^J—m 2 L 2 , the 
dynamical exponent z acquires contributions from two 
more parameters mL and A, namely, z = z(ni{, (3, mL, A). 
Note that for </> = 0, z does not depend on either mL 
or A. For example, without the scalar backreaction, z = 
2 for m,f = 0.36 and j3 — 19.95. Including the scalar 
backreaction, z is modified to approximately 2.025 for 
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m 2 L 2 = —2.2, A = 10 and the same values of mf and /?. 

Since the far-interior solution is of the Lifshitz type, 
following essentially the same argument as before, one 
can analyze equation (48) to show that demanding a non- 
trivial regular solution for <fi(r — > oo) requires the scalar 
to sit at a minimum of V((f>). This is exactly the case 
when z acquires contributions from ml and A. The form 
of this non-trivial regular solution in the far-interior re- 
gion is simply given by (26), with z and g given by their 
new backreacted values. 

In order to flow up to the asymptotic region, we per- 
turb away from the far-interior Lifshitz solution (49) by 
taking the ansatz 

f{r) = jrA 1 + hr a + ■■■), (50) 
g(r) = ±{l + Sl r a + ---), (51) 

/i(r) = 1(1 + fur" + •••), (52) 
0(r) = (f) O + br s + ■■■ . (53) 

We demand the perturbation to grow asr->0 and not 
blow up as r — > oo. We also discard the perturbation 
which gives rise to turning on a finite temperature. Sub- 
stituting the above ansatz back into the equations of mo- 
tion, and taking into account the above consideration for 
the type of perturbation, we then find that the expres- 
sions for fi, 0i, f)i and a take the same form as in the 
case with no scalar backreaction, although their numeri- 
cal values are changed due to the change in z. Similarly, 
the expression for S is the same as 5- in (27). Thus, 
similar to the case without backreaction, the scalar field 
sits at a minimum of V((/>), which results in z acquir- 
ing a contribution from mL and A. The perturbation 
described above can then be thought of as an irrelevant 
deformation of the (2+l)-dimcnsional IR Lifshitz fixed- 
point, which is assumed to be the holographic dual of 
the bulk solution (49) with <p siting at the bottom of 
the potential. 

Note that the backreacted geometry is asymptotically 
AdS4, and the leading-order expansion of <fi (expanded 
around <p = 0) in that region is the same as in (29). Thus, 
if the non-trivial regular solution for (f> in the far-interior 
region could be matched to a normalizable solution in the 
asymptotic region, then the boundary theory dual to this 
background would be in a condensed phase where the 
dynamical exponent z is different from the value when 
the boundary theory is in an uncondensed phase. So, 
one can by tuning a knob, which, for simplicity, could be 
taken to be the conformal dimension of the operator $, 
drive the boundary theory from an uncondensed phase, 
(<&) = 0, with z — z\ to a condensed phase, ($) ^ 0, 
where z — z-i ^ z±. 

In order to find a normalizable solution for (f> in the 
asymptotic region, we scan different values of b (or, more 

precisely, the scale-invariant combination b/f^ a ) to de- 
termine whether there is any solution for which A = 
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FIG. 6. Plot of B versus A for A = 1000. Here z = 4, m f = 0.36 
and m 2 L 2 = —2.2. The BF bound of the far-interior region is 
m 2 c L 2 ~ -1.952. 

while B ^ 0. Figure 6 shows a plot of B against A for 
the case of z = 4 and A = 1000. As the plot shows, 
the qualitative behavior of B against A is similar to the 
case when the scalar backreaction is ignored. The plots 
show that the boundary theory dual to the backreacted 
background is in a phase where ($) ^ 0. 



IV. NEUTRAL SCALAR CONDENSATION AT 
FINITE TEMPERATURE 

In this section we extend our studies of the conden- 
sation of the operator $ to finite temperature. The 
boundary field theory in this case is dual to the finite- 
temperature electron star background [28, 29]. As we 
alluded to earlier, a finite-temperature electron star back- 
ground exists only for temperatures less than a critical 
temperature %. For T > T*, the background is a non- 
extremal RN-AdS4. The condensation of a neutral scalar 
operator for the theory dual to a non-extremal RN-AdS4 
background was already addressed in [8] . Hence, the tem- 
peratures we consider in this section are always in the 
range T < T* for which the electron-star background is 
the preferred solution. 

The operator <j> is dual to the neutral scalar field (j) 
whose bulk action is given in (19). The assumption re- 
garding the mass m of the scalar field is as before. Also, 
we first consider the case where the scalar backreaction 
is ignored, namely, for large A. The equation of motion 
for the scalar field is given in (22). We are after a non- 
trivial solution which is regular as r — > r$ and normaliz- 
able as r — > 0. Such a solution in the bulk characterizes 
a phase of the boundary theory where the dual opera- 
tor $ condenses. At finite temperature, demanding a 
non-trivial regular solution for the scalar <f> close to the 
horizon does not require the scalar to be sitting at the 
bottom of the potential. This can easily be verified by 
expanding equation (22) close to the horizon. Note that 
close to the horizon, in the so-called inner region of the 
finite-temperature electron star background, the metric 
functions are given by (17). Indeed, once expanded close 
to the horizon, equation (22) admits a regular solution of 
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FIG. 7. Plot of B versus A for m 2 L 2 = -2.2 and T/T c ~ 0.443. 
For the background, we set mf = 0.7 and /3 = 10. 
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FIG. 8. T c versus m?L 2 . For the background, we set mf = 0.7 
and /3 = 10. The shaded region represent the condensed phase, 
and the dashed line marks the BF bound of the far-interior region 
as T -> 0. 

the form 

4>{r) = 0o + M r a -r) + 4> 2 {r - r Q ) 2 H , 

where <fri, <f>2 and all other coefficients are determined in 
terms of (f) - 

The asymptotic behavior of </> (expanded around 4> — 
0) will be given by the same expression as in (29), 
as both the zero- and finite-temperature electron star 
backgrounds are asymptotically AdS 4 . As in the zero- 
temperature case in the previous section, we choose to 
work in the standard quantization where A is the source 
and B is the vev of the dual operator $. Thus, by 
varying <p , one can search for a normalizable solution 
asymptotically with A = and B =/= 0. For numerical 
computations in this section, we define a dimensionless 
radial coordinate by rescaling r — > r$r, and also work 
with dimensionless coefficients A and B defined through 
A — > r^~ 3 A and B — > r^ A B. Our numerics show 
that the operator $ condenses below some critical tem- 
perature T c as long as the mass squared of the scalar 
field satisfies -9/4 < m 2 L 2 < m 2 c L 2 . Note that the 
lower limit is there to ensure stability in the asymp- 
totic AdS 4 region. Figure 7 shows a plot of B versus 
A for m 2 L 2 = —2.2 on an electron star background with 
m f = 0.7 and f3 = 10. For the plot, T/T c ~ 0.443. No- 
tice that as 0o — > ~4>0i A ~A and B — > — B. So, 



without loss of generality, we have just plotted half the 
curve of B versus A in Figure 7, the half for which (f>Q 
is non-negative. Our numerics also show that T c — > as 
m 2 L 2 approaches the critical value from below, as shown 
in Figure 8, although, due to numerical difficulties, we 
could not probe larger values of m 2 L 2 , namely the ones 
very close to m 2 L 2 . Moreover, as shown in the plots of 
Figure 9, we find that close to T c , B ~ (1 - T/T c f" with 
/3 C ~ 0.506, dB/dA\ A=0 ~ (1-T/T c )-t<= with 7c = 1.015, 
and at T = T c , B - A 1 /" 5 - with S c ~ 3.07. Also, Fig- 
ure 9(d) shows a plot of the change in the free energy 
density AO = ^condensed - O uncondonsc d versus the tem- 
perature. It shows that whenever there is a condensed 
phase, its free energy is always lowered when compared 
to the uncondensed phase. Hence, the condensed phase 
is energetically favored over the uncondensed phase. As 
temperature approaches T c from below, the difference in 
free energy approaches zero as AO ~ (1 — T/T c ) Vc , with 
the numerically obtained value of v c ~ 2.061. 

The numerically obtained critical exponents /3 C , j c , S c 
and v c are indicative of a second-order phase transition 
with mean field exponents. The finite-temperature phase 
transition obtained here is indeed similar to the result 
found in [8] for the scalar condensation in the boundary 
theory dual to the non-extremal RN-AdS4 background. 
This is perhaps not surprising as the nature of this phase 
transition at finite temperature could be traced back to 
the analyticity of A and B as a function of 4>o close to 
the horizon [8]. 



A. Backreaction at Finite Temperature 



We now take into account the backreaction of the neu- 
tral scalar field on the finite-temperature electron star 
background. Again, since the neutral scalar does not 
couple directly to the fermion, the fermion density of 
states is given by (7) . This assumes that the effect of the 
temperature on the fermionic density of states is negli- 
gible. We also assume that the fermion dynamics are 
determined by the local chemical potential, with a form 
given in (8), with the same formulae for a, p, p given as 
in the case without scalar backreaction. Similar to the 
case without scalar backreaction, the finite-temperature 
background is divided into three regions. The radii ri 
and r 2 which separate the three different regions are so- 
lutions to the equation fi\ oc (r) = m f . Our discussion 
of the backreacted background in the following will be 
brief as its characteristics are very much similar to the 
case without the backreaction of the neutral scalar. In 
the inner region, r > r > r 2 , the solution is obtained 
by solving the equations (44), (46), (47) and (48) with 
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A T/T, 

FIG. 9. Plots of (a) B versus T/T c , (b) dB/dA at A = as function of 1 - T/T c (c) B versus A at T = T c and (d) Af2 versus T/T c . We 
set m 2 L 2 = — 2.2 and for the background, m = 0.7 and /3 = 10. From these plots, the critical exponents /3 C , 7c, <5 C , and u c are found to 
be fi c ~ 0.506, 7c — 1.015, <5 C ~ 3.07 and t> c ~ 2.061 which indicate a second-order phase transition with mean field exponents. 



cr = 0, namely 



1 ff(r) | g'(r) 
r V f( r ) 9{r) 
f'(r) h'{rf 



2A 



r 2 4A v 



2/(r) 

2 



+ 



h"(r) + ^rft'(r)^(r) 2 , 





= 0"(r) + 



1 //'(r) S'(r) 



2 V /W ff(r) 



- 5 (r)L 2 V"(0). 



(54) 

(55) 
(56) 

(57) 



Our objective in this subsection is to study the conden- 
sation of the operator <f> in the boundary theory dual 
to the backreactcd finite-temperature electron-star back- 
ground. This means that we are after a solution for <p 
which is regular at the horizon (as well as normalizable 
in the asymptotic region). So, we start by taking the fol- 
lowing ansatz for the expansion of the solution near the 
horizon 



9(r) = 



9o 



+ 9i + ■ ■ ■ 

Mr 



r-r 
f(r) = fa(r-r a ) 
h(r) = h (r - r ) + h^r - 
4>(r) = 0o + 4>i(r - r ) + ^{r - r ) 2 + ■ 



-r f + - 



(58) 

(59) 
(60) 
(61) 



It is more convenient to rescale the radial coordinate such 
that the horizon is at tq = 1. Substituting the ansatz 



back into the equations (54)-(57), one finds that all of 
the coefficients in the above expansion are determined in 
terms of ho, fo and 4>o- Note that the value one chooses 
for fo fixes the normalization of time. In the intermedi- 
ate region r 2 > r > r x where the fluid energy density, 
pressure and charge density are all non-zero, one has to 
solve the equations (44), (46), (47) and (48). In the outer 
region r < ri, the space-time is given, once again, by 
the solution to the equations (54)- (57). The outer-region 
solution is asymptotically AdS4, with the leading expan- 
sion of the scalar field given by (29). Note that in the 
backreacted case, the temperature of dual field theory 
is given by (47rc)~ 1 \df (r )/dr y/l/Jjrpjg (r ) | , where c is 

given by \/r 4 f(r)g(r)\ r ^o- Also, the chemical potential 
of the dual theory is given by h(r -4- 0) /c. 



The condensed phase (in the standard quantization) 
is once again characterized by A = and B ^ 0. The 
plots in Figure 10 show the behavior of B versus A for 
different values of the coupling A. As the plots show, the 
condensed phase persists even as A decreases, where the 
backreaction of the neutral scalar field becomes more im- 
portant. As a test of our numerics for A — > oo, the curves 
produced with the scalar backreaction present approach 
the curve obtained for the case without the scalar backre- 
action. We have checked that even if the scalar backreac- 
tion is included, the finite-temperature phase transition 
between the condensed and uncondensed phases is still 
second-order with mean-field exponents. 
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FIG. 10. Plots of B as function of A. The black curve corresponds 
to the case without the backreaction of the scalar field, while the 
blue, green and red curves correspond to the scalar backreaction 
with A = 10, 20 and 100 respectively. The mass square of the scalar 
is m 2 L 2 = —2.2 for all the curves. For the non-backreacted case, 
T/,u ~ 0.00045 with m f = 0.7 and /3 = 10. 



V. DISCUSSION 

The Hertz-Millis [3, 4] theory is an attempt to describe 
quantum criticality in a metallic system in the presence of 
some type of magnetic order described by a neutral order 
parameter. In obtaining the low-energy theory in terms 
of the bosonic neutral order parameter, the fermions, 
which are gapless, are integrated out. This procedure 
has no guarantee of working since the fermions belong to 
the low-energy sector. In d = 2 + 1, this procedure is fa- 
tal as an infinite [5-7] number of marginal operators are 
generated by the integration procedure. Holography ob- 
viates integrating out the fermions because they can be 
included non-perturbatively in the background. While 
the geometry remains AdS4 asymptotically, the fermions 
deform the space-time in the interior giving rise to an IR 
Lifshitz fixed-point in which the dynamical exponent is 
inherently finite. What we have shown is that the con- 
densation of the neutral scalar operator in the bound- 
ary theory dual to the electron star background is con- 
trolled by the IR Lifshitz fixed-point whenever the mass 
square of the neutral scalar field in the bulk violates the 
BF bound. Figure 1 shows explicitly that the boundary 
theory undergoes a quantum phase transition simply by 
tuning the mass square of the bulk scalar field. It is the 
condensation of the scalar that gives rise to the new IR 
scale, see equation (43), that depends exponentially on 
the distance from the critical point. Consequently, the 
underlying T = transition is of the BKT form. At 
finite temperature, we find that the phase transition is 
second-order and is described by mean-field exponents. 
Backreaction of the condensed scalar on the geometry 
leads to a change in the dynamical exponent across the 
phase transition. This is a particularly attractive feature 
of this theory as there are a variety of systems [31, 33] 
in which the condensation of a neutral field changes the 
dynamical exponent across the phase transition. 

There are a number of immediate applications of our 
work. First, our model can be used to described anti- 



ferromagnetic phases in condensed matter systems with 
finite dynamical critical exponents. To do so, we would 
simply follow the analysis discussed in [8] by embedding 
the neutral scalar field into a triplet charged under an 
SU(2) gauge symmetry in the bulk; see also the discus- 
sion in [24]. This bulk SU(2) symmetry are to model the 
SU(2) spin symmetry in the dual boundary theory. The 
antiferromagnetic ordering then corresponds to the spon- 
taneous breaking of the SU(2) spin symmetry to U(l) in 
the low energy limit. To holographically model the tran- 
sition to the antiferromagnetic phase, where there is a 
staggered spin order parameter with zero spin density, 
we introduce in the bulk an SU(2) gauge field A a M along 
with a triplet </> a charged under the SU(2) gauge group. 
Here, a = 1,2,3. (Note that 4> a is neutral under the 
U(l) gauge group.) A a M is dual to the spin density in 
the boundary theory which will be set to zero, while 4> a 
is dual to the staggered order parameter. By embedding 
the neutral scalar field <f> into </> a , we can explicitly break 
the SU(2) spin symmetry when the mass of the scalar 
field falls below the BF bound of the far-interior Lifshitz 
region. Further, by studying perturbations around the 
symmetry breaking solution, one can find that there are 
two gapless Goldstone mode with linear dispersion rela- 
tions. The argument presented in [8] is independent of 
the form of metric and can be applied to electron star 
solution also. 

Second, we could also model nematic order. As men- 
tioned in the introduction, nematic order [2] is an in- 
stance of a neutral tensor condensate. The order pa- 
rameter for a nematic phase is a traceless symmetric 
second-rank tensor. Indeed, one way for holographically 
constructing such a phase would be similar to the holo- 
graphic setups studied in [40-42] for d-wave supercon- 
ductivity, except that the order parameter here is neu- 
tral. Thus, one introduces a massive neutral traceless 
symmetric rank-two tensor field in the bulk. A rotation 
in the (x, y)-plane can be used to set one component of 
the tensor field equal to zero. The equation of motion 
for the remaining component, after a field redefinition, 
will be identical to the equation of motion for a neutral 
scalar field. So, the analysis for the condensation of the 
dual tensor operator in the boundary theory is essentially 
mapped to the analysis for the condensation of a neutral 
scalar operator that we presented in this paper. It would 
be interesting to see the Pomeranchuk instability [43, 44] 
in this holographic setup as the quantum critical point is 
approached. Of course, this would require coupling the 
tensor field to a fermion. One should look at our heuristic 
description of a holographic quadrupolar nematic phase 
from an effective field theory point of view, mainly be- 
cause, as of now, there is no consistent action describing 
the dynamics of a massive spin-two field which would 
be causal and free of ghost in a space-time which is not 
Einstein; see [42, 45] for more discussions. 

Third, there are numerous examples in condensed mat- 
ter as well as particle physics (most notably in QCD) in 
which several ordering tendencies compete. In the in- 
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stances in which competing order involves neutral con- 
densates in metallic systems, the holographic method 
used here can be tailored to capture the relevant physics. 
The Z 2 symmetry in the case of a single neutral scalar 
field can easily be enlarged by including more scalar fields 
in the action (19) and then choosing the interactions to 
be invariant under the desired symmetry. Depending on 
which dual operators condense, a variety of interesting 
symmetry breaking patterns could be potentially realized 
in the boundary theory. 

Finally, another extension would be the inclusion of 
probe fermions coupled to the neutral scalar field. For 
the RN-AdS4 background, fermion correlators in the dual 
boundary theory have been studied extensively; see, for 
example, [24, 41, 46-60] . A similar kind of analysis can be 
done here for the electron star background. In fact, the 
results on the electron star background would in principle 
be more relevant to condensed matter systems because 
the IR theory is controlled by a Lifshitz fixed-point and 
moreover, the entropy density of the system goes to zero 
at zero temperature. While some work along these lines 
has been initiated in [61], probe fermions which couple to 
an order parameter, either charged or neutral, have not 
been included. Including such couplings could give rise to 



a non-trivial dispersion for the fermions around the Fermi 
surface. One interesting effect that could be studied in 
this approach is whether or not a neutral condensate with 
a non-trivial momentum structure could, upon coupling 
to the probe fermions, give rise to a Fermi arc structure 
indicative of the pseudogap phase of the cuprates. This 
is particularly appealing since most models [62] of the 
pseudogap involve some type of neutral order parame- 
ter. Recent holographic constructions [41, 57] of Fermi 
arcs are not particularly relevant to the pseudogap phase 
of the cuprates because in these construction the exis- 
tence of the Fermi arcs are due to the coupling between 
the probe fermion and a charged superconducting order 
parameter. 
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